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Abstract 

Wc present an original study on the numerical stabiliy of explicit schemes solv- 
ing the incompressible Euler equations on an open domain with slipping boundary 
conditions. Relying on the skewness property of the non-linear term, we demonstrate 
that some explicit schemes are numerically stable for small perturbations under the 
condition St < C6x^^^^^^~-^^ where r is an integer, 6t the time step and Sx the space 
step. 

1 Introduction 

In order to achieve stability of a numerical scheme solving the incompressible Euler equa- 
tions in a divergence-free dicretisation frame, we seek for a criterion that link the time 
and the space steps. During unsteady incompressible fluid simulations with the help of 
divergence-free wavelets [2], we observed the CFL-like condition: 6t < C5x^^^ for a cen- 
tered scheme order two (4.2), presented in [3, 2]. Either bibliographical data don't give a 
satisfying explaination to this phenomena [5, 4], either it relies on a Von Neumann sta- 
bility analysis for the advection equation where velocity is a constant vector [7], which is 
not our case. 

This work takes advantage of ideas presented in R. Temam [5, 4] for its context: 
incompressible Euler equations with the use of the skewness property (lemma 3.1). But, 
at the difference with references [5, 4], we are not in the context of finite elements, but 
assume a divergence-free space of discretisation. Meanwhile, we make the connection with 
the results of Von Neumann stability for the convection equation [7] . 

In the following Note, we manage to derive the propagation in norm of a small 
perturbation e„ of a regular solution u„ of the numerical explicit scheme 2.1 to solve the 
incompressible Euler equations (1.1) in a relatively accurate way. 

In this case, we establish a stability criterion of the form 5t < C5x'^^^^'^'^~^^ with r an 
integer and C an explicitly computable constant, for various explicit schemes. Contrarily 
to Von Neumann stability [7] which assumes infinite or periodic domain and constant 
advection velocity, we establish this stability condition solely under regularity assumptions 
on the velocity. 
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The Euler equations modelise incompressible fluid flows with no viscous term: 

-^ + (u- V)u- Vp = 0, divu = (1.1) 
at 

The use of the Leray projector P which is the L^-orthogonal projector on the divergence- 
free space, aflows us to remove the presure term: 

^ + P[(u.V)u]=0 (1.2) 

We start from a discretisation u„ of the solution u in time and in space. Then we consider 
a perturbed solution u„ + e^. Assuming the regularity of u and the consistency of the 
scheme, we are interested in the evolution of the perturbation e„ for different explicit 
schemes in time. 

2 Discretisation in time and in space 

In order to solve equation (1.2) numerically, we discretise the solution u„ = u(nSt) in a 
divergence- free space Vdivo{6x). Such spaces appear in spectral codes using the Fourier 
transform [1], or are produced in a stable way thanks to divergence- free wavelets [6, 2]. 
The parameter 6x stands for the smallest space step in the discretisation space Vdivo('^2;). 
Hence, every function u„ in Vdivo{Sx) satisfles: 

||SjU„||LP < C{p)5x~'^\\Un\\LP 

where 1 < p < -|-oo, and 5jU„ denotes a partial space derivative of u„,. The constant C{p) 
can often be taken equal to 1, and we will do so in the following. 

Concerning the discretisation in time, we consider different explicit schemes proceeding 
in several steps, like Runge-Kutta schemes: 

e-i e-i 

U„(0) = Un, U„(f) = ^ an U„(i)-X] [(""(0 ■ '^)"n(j)] iov I < £ < k, U„+i = U„(fc) 

i=0 i=0 

(2.1) 

where P stands for the orthogonal projector on the discretisation space Vdivoi^^)- One 
can notice that PoP = PoP = P. 

3 L^-stability condition for a small perturbation 

Actually, the small error e„ that we introduce corresponds to oscillations at the smallest 
scale in space Vdivo{Sx). This stability error propagates and may increase at each time 
step. In what follows, we demonstrate that under some precise CFL-like conditions, the 
norm of this small error e„ is amplifled such that: 

||e„+i||L2 < (1 + C<5t)||e„||i2 (3.1) 

where C is a constant that neither depends on 6x nor on 6t. 

Thus, after a time elapse T, the error increases at most exponentially as a function of 
the time: 

\\et,+T\\L^ < {I + C6tf/^' \\et,\y < e^'^Wetoh^ (3.2) 

For the following stability study, we will need the skewness property of the trans- 
port term. This property is utilized for the stability of the incompressible Navier-Stokes 
equations in [4]. 
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Lemma 3.1 Let u, v, w € H^{Q,)'^, H^{Q,) denoting the Sobolev space on the open set 
n C R'^, be such that (u • V)v, (u • V)w G L^. //u G H^wfliQ) = {f G {L'^{n)))'^, div f = 
0}, then 

< V, (u • V)w >L2(n)= - < (u • V)v, w >i2(n) 
Corollary 3.1 With the same assumptions as in lemma 3.1, 



< V, (u • V)v >L2(n)= / V • (u • V)v(ix 





In the scheme (2.1), we denote by £^(1) the stability error at level I. Then, under 
the condition 5t = o{5x) and for e„ small enough, most of the terms appearing in the 
expression of become negligeable compared with: 

• the terms 6t^ F^i^n) where F{en) = IP[(u„ • V)en] and F^ = FoFo---oF,i times. 

• the term 6tF[{en ■ V)uJ, 

Hence, for the scheme (2.1), we find the following expression for e„(£): 

e 

en(i) =Y,^e^SfF\en) + PeSt¥[{en ■ V)u„] + o() (3.3) 

where function o() gathers the negligeable terms. 

i-i i-i i-i e-i 

an = ^ aejOji + ^ bijOji^i and Pe = '^ aij(3j + ^ bij 
j=i i=i-i i=i j=o 

As a result, we find the following expression for e„+i: 

k 

En+i = ai6fF\en) - 6tn{en ■ V)uJ + o() (3.4) 

Starting from this expression and using the fact that according to lemma 3.1, 

if i + j = 2£ + 1 for ^ G N 

Il2 



<F\er,),F^ien) >LHn)= { 



(_l)^-*||i?^(e„)||2 if i + j = 2i for ^gN 



(3.5) 

we compute the norm of £n+i as a function of the norm of e^: 

k 

\\en+i + Stn{en-V)un]\\l2=Y.Se6t^'\\F\en)\\l2+o{) (3.6) 

e=Q 

with 

min{e,k-e) 

St= ^ {-ly ae_jai+j (3.7) 

j=-mm{£,k-e) 

For consistency needs of the numerical scheme, we must have 5o = 1. If, on an other hand 
we suppose Si = S2 = ■ ■ ■ = Sr-i = and 5^ > 0, knowing that in the discretised space 
Kiivo('^2;), 

\\F"'{£n)\\L2 < ||un||£oc-^i^^|^ and ||P[(e„ • V)u„]||j^2 < ||e„||i2 ||Vu„||i,oo 



3 



we derive: 

( ( 5t^'''~^Sr \ \ 

||en+l||L2 - ( ^ + ( l|VUn||L-= + 2fa^^ ^^^'^ ^°^^) "^""-f^^ ^^'^^ 

If, on an other hand we assume the consistency, there exist constants Aq and Ai such that 
llunllioo < Aq and ||Vu„||ioo < Ai when 5x and 5t go to 0. Hence, the numerical scheme 
(2.1) is stable for small pertubations under the condition: 

5t < C6x^ (3.9) 
That brings the following theorem out: 

Theorem 3.1 An order 2p scheme solving the incompressible Euler equations is numer- 
ically stable for small perturbations at worst under the CFL-like condition: 

6t < C6x^ (3.10) 
Proof: For an order 2p scheme, we have the following equality, point by point: 

U„+i = U„ + 6tdtUn + ^dttUn H h J^^^t^^n + o{6t^P) 

Considering dfUn = F[(u„ • V)u„] and introducing a small perturbation leads to: 

Sn+i = en+dtFien) + —FoF{en) + - • ■ + -—FoF^P~\en)+6tF[{en-V)un]+o{) (3.11) 

2 {2py. 

with F{e) = P[(u„ • V)e] and o() gathering the terms that are negligeable under the 
condition 6t = o{6x). Then for q G [l,p]. 

Which allows us to conclude, as stated at line (3.8). 



Remark 3.1 ^ Von Neumann stability analysis would have proceeded as follows. We 
compute the evolution of the Fourier mode ip{n6t) = ipne^'^'^ with ( G M'^, for the advection 
equation dtip = —a • Vtp, with a = u a constant velocity. As Vip = iC, ■ ip, for the scheme 
(2.1), we find v3„,+i = S,ipn with, as for computation (3.4), 

k k 
^ = Y^aj{-ia-Cy and then \^\^ = Sf^^l^ ■ C\^^ 

j=o e=o 

The coefficients Si have the same expression ( 3. 7) as when we used the skewness property. 
As, on an other hand, we have |a • C| < ||a||/(5x in the discretisation space Vdivo('52;), we 
find the same stability criterion if we want to have |^j < 1 + C5t. 
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4 Examples 

Let ^0 = sup4g[o,T], xen |u(i,x)| and Ai = sup^gp ^.^^ |Vu(t,x)|. We propose to apply 
our stability analysis to some classical schemes. 

The simplest example is the Euler explicit scheme, order one in time: 

u„+i = u„, -6tF [(u„ • V)u„] (4.1) 

For this, we find: 

\\e^^,\\^,<h+{^— + Ai)5t\\\en\\L^, and the CFL : 5t<2Ci — ] 

An improved version of this scheme allows us to construct an order two centered scheme: 
I U„+i/2 = u„ - f P [(u„, • V)u„] 

[ U„+i = Un-dt¥ [(u„+i/2 • V)u„+i/2] 

For this scheme, stability is slightly improved: 

||e„+i|Ii2 < (^1 + + IknllL^ hence the CFL : 5t < 2C^/^ (^^^ ' 

For Runge-Kutta scheme of order 4: 
f u„(i) = u„ - f P [(u„ • V)u„] 

U,i(2) = U„ - f P [(U„(i) • V)u„(i)] 
U,i(3) = U„ - (5t P [(U„(2) • V)U„(2)] 
I U„+i = U„ - f P [(U„ • V)U„] - f P [(U„(i) • V)u„(i)] - f P [(U„(2) • V)U„(2)] - f P [(U„(3) • V)u„(3)] 

(4.3) 

theorem 3.1 allows us to predict a CFL-like condition bt < C5x^^^ at worst. Computations 
show that actually: 

5i = 52 = and 53 = - — , ^4 = — 

72 bib 

Hence our study doesn't fully apply to this case. 

The order two Adams-Bashford scheme doesn't remain to the definition (2.1). It goes 
as follows: 

u„+i = Un- -5t P [(u„ • V)u„] + -5t P [(u„_i • V)u„_i] (4.4) 

Nevertheless, it is possible to perform computations similar to those of part 3. First, as it is 
order two, according to theorem 3.1, it is stable at worst under the condition 5t < C6x^/^. 
Further computations show that: 

/ ^^4 \ / 5x\ 

lkn+i||L2 ^ ( 1 + Js^^o + ) 11^" 11^' inducing the CFL: 6t < 2^1^ f — j 

As a conclusion, many usual schemes for simulating the fluid flows verify a stability 
condition of the type bt < C6x'^'^^^'^^~^^ with r an integer. 

On can remark that if we have relation (3.11) at order m with no term of the type 
F o ■ ■ ■ o F{en) in o() (which is, for instance, what Runge-Kutta schemes satisfy), then the 
related scheme will have to verify a CFL-like condition of the type 6t < C'fe^™'^^)/'" if 
m = 1[4] and 6t < if m = 2[4]. 
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